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Abstract. The expected range of a sample of n + 1 normally distributed 
variables is known to be related to the mean width of a regular n-simplex. We 
show that the expected maximum fj, n of a sample of n half-normally distributed 
variables is related to the mean width of a regular n-crosspolytope. Both of 
these relations have mean square counterparts. An expression for fx 5 is found 
and is believed to be new. 

Let C be a convex body in R™. A width is the distance between a pair of parallel 
C-supporting planes (linear varieties of dimension n — 1). Every unit vector u G W 1 
determines a unique such pair of planes orthogonal to u and hence a width w(u). 
Let u be uniformly distributed on the unit sphere S" 1-1 C M n . Then w is a random 
variable and 

E (ws) = — arccos I - 
7T \3 

for C = the regular 3-crosspolytope (octahedron) in M 3 with edges of unit length. Our 
contribution is to extend the preceding mean width result to regular n-crosspolytopes 
in M. n for n < 6, starting with an integral formula in [Tj [2]. We similarly extend the 
following mean square width result: 

E (^ 3 2 ) = 

which, as far as is known, first appeared in [3]. 

1. Order Statistics 
Let X 1 , X 2 , X n denote a random sample obtained by "folding" Normal (0,1) 
across the origin, that is, with density function / and cumulative distribution F: 

m = \jl exp (-y) , F{x) = J /(£) d£ = erf ^ 

o 

for x > only. 
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The first two moments of the half-normal maximum 

m n = max{Xi, X 2 , X n } 

are given by [I] 

oo oo 

H n = E(m n ) = n J x F(x) n ~ 1 f(x)dx, v n = E(ra£) = n J x 2 F{x) n ^ 1 f{x)dx. 
o o 

For small n, exact expressions are possible [U E] : 



^ = ^ = 1.128.., 

^3 = ^ = 1.326..., 

M4 = ^(1 - 4£ a ) = 1.464.., 

^ = 5 (-V2 + 8S 2 + 16\/2T() = 1.569... 

/i 6 = || (1 - 8S 2 + 16T 2 ) = 1.653... 



^2 


= 1 + 


2 = 1.636.., 

TT ' 


"3 


= 1 + 


^ = 2.102.., 

TT 

= 2.470.., 


l/ 4 


= 1 + 




= 1 + 


20 ^ (1 453) - 2.773 



where [7] 



5V 



tt/4 



dx 



1 



- 7 v/&T 

v 



V2 



tt/4 tt/4: 




dx dy 



sec(x) 2 27r 

TT S2 



arcsec (k + 1) 



J J J2 + sec(x) 2 + sec(y) 2 27r 
v 



arcsec I 1 + 



6 



2 — tan(-z) 



7r/4 tt/4 




dx dy 



ir* j j 1 + sec(x) 




sec(y) 1 



1 - 2V2S 2 1 7 2 arCS6C ( 2 ~ 



16 



2y/2n 2 



l+4cot(,z) 2 



1 - 



6 



dz. 



2-tan(z) 2 



The preceding table complements an analogous table in [8] for first and second mo- 
ments of the normal maximum (rather than half-normal). Similar expressions for 
v% = 3.032... and fi 7 = 1.723... remain to be found. Our expression for /x 5 is believed 
to be new and details will be provided shortly. 
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2. Starting Point 
A regular n-crosspolytope with edges of unit length has mean width [TJ [2] 

oo 

\ 2 J Q 

Integrating by parts yields 



oo oo 



e~ 2x2 erf (x) n ~ 2 dx = / x e~* eifix)^ 1 dx 

n-lj 

o o 

OO j 

= , — r / xe _:r2/2 erf ( ) dx 
o 

oo 

= - * jxf{x)F{x) n - l dx. 
2y2\n — 1)J 
o 

We recognize the latter integral as fi n /n; hence 

F (rt) 

and conjecturally 

E (w 2 n ) = -u n . 
n 

Numerical confirmation for n < 6 is possible via the computer algebra technique 
described in [3]. 

In summary, we have mean width results 

4 

E(w 2 ) = - = 1.273239544735162..., 

7T 

E(iu 3 ) = 6S 2 = 1.175479656091821..., 
16 

E (wA = — (1 - 4S 2 ) = 1.101845693159859..., 

7T 

15 / \ 
E(w 5 ) = — (^-^+8^2 + 16V2T^J = 1.043421681509785..., 

32 

E (w 6 ) = — (1 - 8S 2 + 16T 2 ) = 0.995378656038812... 

7T 
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and mean square width results 

E {wfj = 1 + - = 1.636619772367581..., 

E (wl) = - | 1 + — ) = 1.401771860562389..., 

3 \ 7T / 

E (wl) 4 1+ ^ I = 1.235105193895722..., 
2 \ on I 

E (wj?) = 1 (\ + (l _ 453) \ = 1.109499626837713.... 

3. Details Underlying /i 6 

To compute n Q is less difficult than to compute /i 5 . One approach is to use [9l[l0j[TT], 
relating moments of order statistics for a symmetric distribution with moments of 
the corresponding folded distribution. We will follow [7], however, because no other 
approach seems to be applicable to evaluating fi 5 as well. 
As a prelude, 

</« = ->) 

hence 

00 00 00 

/i 6 = 6 J tF{tff{t)dt = -6 J F{tfj t f{t)dt = 30 J F{tff{tfdt 
000 
via integration by parts. Also 

F(t) 2 = —J J exp ( ^ - j dxdy = — J J exp ^— — ^ r dr d6 

00 00 



via transformation to polar coordinates. 
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We deduce that 



30 



oo r tt/4 

/Mm- 

o L o 

oo oo vr/4 

//w* -|//« 



t 2 sec(ey 



d9 



t 2 sec(6f 



f(t) 2 dt 



f(t) 2 dddt 



+- 



16 





oo tt/4 tt/4 



71" 




exp 



t 2 (sec(y?) 2 + sec(-0) 5 



f{t) 2 dpdipdt 







oo tt/4 tt/4 

1 8 32 f f f 

^~7^ 2+ w/ /exp 





7r/4 7r/4 



32 hi 



a/tt a/tt tt 3 V 2 J J a/2 + sec(y2) 2 + sec(^) 



t 2 (2 + sec(y?) 2 + sec(» 2 ) 



dip dip 



dip dip dt 



o o 



1 8 32 [rTn 2 1 

— /= 7=^2 H — a W 7T^ T 2 = — f= (1 - 86 2 + 16T 2 ) . 

y/n 7T 3 V 2 a/2 a/tt 



As an interlude, 



tt/4 

/_ 

J J~K 



dip 



^Jk + sec(» 2 2 \/« 



7T 1 

S K = - p= arcsec + 1) 



hence 



7r/4 7r/4 



tt/4 







a/2 + sec(y?) 2 + sec{ip) 2 2 J \/ 2 + sec(y?) 



. arcsec (3 + sec(y?) 2 ) c?<^ 



upon substituting 2 + sec (79) 2 for /t. Now let 



then 



which implies that 



2sin(v?) 2 = 3sin(#) 2 , 
2-2cos(v?) 2 = 3-3cos(^) 2 



sin(#) 2 = - sin(v9) 2 , cos(#) 2 = - (1 + 2 cos(y9) 2 ) 
3 3 
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hence 

4sin(</?) cos(ip)dip = 6 sin(#) cos(9)d9 

= -a/(2 sin(^) 2 ) (1 + 2 cos(v?) 2 ) d9 
= 2 V2 sin(v?) \A + 2 cos(<^) 2 d0 

hence 

dp = i^?^ = -^ v / 2 + sec(^ ^. 
v/2 cos(v9) v^ v 

The preceding integral becomes 



since 

2 + sec(v?)2 = 2 , 1 _3-2sin(^_3-3sin(^_ 6cos(^ 



l-sin(v?) 2 l-sin(v?) 2 l-fsin(^) 2 2-3sin(0) 2 
6cos(#) 2 6 



2cos(#) 2 -sin(#) 2 2-tan(#) 2 ' 
Finally, as a postlude, 

2arccos(C) = arccos (2( 2 — l) 
and if = 7r/4 implies that 2 = 3cos(6>) 2 , thus the upper integration limit is 



9 = arccos ^\Z§^ = I arcc °s (| — l) —\ arcsec (3) = n S 2 . 

4. Details Underlying // 5 
Integration by parts is not employed this time. Beginning with 

oo 

^ = ^JtF{tff{t)dt 
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it follows that 



5 



oo r tt/4 

f , 4 f ( t 2 sec(0H 
o L o 

oo oo tt/4 

J t f{t)dt-^J J exp 



-■ 2 



t 2 sec(#) 2 



tf(t)dt 



tf(t)dOdt 



+- 



16 



o o 

OO 7T/4 7r/4 




cxp 



t 2 (sec(<^) 2 + sec(^) 2 ) 



ooo 



2 2^2 -%S 2 16 [2 
T + 

77 V7T 



oo tt/4 tt/4 



7T 2 V 7T 



i exp 



\= { V2 - (2^2 - 8S 2 



000 

7r/4 7r/4 



16^2 



77 z 



tf(t)d<pdipdt 



t 2 (1 + sec(y?) 2 + sec(^) 2 ) 



dip dip 



dip dip dt 



As an interlude, 



y= (-V2 + 8S2 + 16V2T?) 







1 + sec(v9) 2 + sec(ip) 2 



tt/4 



dip 



7T arcsec (1 



k + sec(ip) 2 An 2k\/k + 1 



hence 

7t/4 7r/4 



tt/4 







1 + sec(y?) 2 + sec(^) J 



77 



arcsec ( 1 



2 



t/4 /, 

_2 /• arcsec 1 + T 

£-(l-2V2S 2 )- / * 

16 J 2(l + secM 2 ) v 



4(1 + sec(v?) 2 ) 2 (1 + sec(y?) 2 ) ^2 + sec(y?) 2 

■.dip 



dip 



-sec(<^) 2 



(1 + sec(y?) 2 ) ^2 + sec(^) 1 



upon substituting 1 + sec((/?) 2 for k. Now, as before, let 

2sin(v?) 2 = 3sin(#) 2 
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then the preceding integral becomes 



2^2 



o 

since 



^Si arcsec ( 1 + - 2 g ) arcsec ( 9 

f V -^ir^wj ^ = 1_ /■ arCSeC ^ 

•' ~~ ^ + 2-tan(0)2 -V^ •/ 1 — 2-lanltf)- 



l+4cot(0) 5 



6 



1 2 = ~ + ^^f) + 2 = 3- 2 , ta 6 n(e)2 = 6-3tan(g) 2 -6 

-!+23^ "1 + 23^ + -2 + tan(^ + 6 

-3tan(#) 2 3 
4 + tan(#) 2 ~ ~1 + 4cot(fl) 2 ' 

Everything comes together at this point, completing the derivation. 

5. ASYMPTOTICS 

We turn now to the asymptotic distribution of m n as n — > 00. Define a n to be the 
positive solution of the equation [121 EE] 

2tt a 2 exp (a 2 ) 

that is, 

1 ln(ln(n)) + ln(47r) 



n 2 , 



in terms of the Lambert W function [13] . An adjustment is needed that corresponds 
to folding (truncating). Define [12] 



, 1 , . . 1— — — — 1 ln(ln(n)l + ln(7r) 

a' n = a n + - ln(2) ~ v^hrR - - 1 7== 1 l , 

then the required density is 

lim — P fv/2 \n(n)(m n - a') < y) = exp (-e~ y ) . 
n->oo ay \ J 

A random variable Y, distributed as such, satisfies 

E(y) = T , E(F 2 ) = ^ + 7 2 

o 

where 7 is the Euler-Mascheroni constant [15J. This implies that 

/ 7 / , 7 s 1 ln(ln(n)) + lnfvr) - 27 



Mean Width of a Regular Cross-Polytope 



9 



and hence 



r (-) 




a' + ' 



ln(n) ln(ln(n)) + ln(7r) - 27 
n A^/n ln(ra) 

More terms in the asymptotic expansion are possible. The first term of this approxi- 
mation is consistent with |16j . 

If we rescale length so that the inradius of the n-crosspolytope is 1 and denote 
adjusted width by w n , then 



E(w n ) ~ V2^ ■ rsj J 2 ln(n) 

V n 



because, in our original scaling, the inradius is y/l/(2n). 

After completing this work, we discovered an article [17] on metric properties of 
the only three regular polytopes in 5-dimensional space or higher. Our paper and 
its predecessor [3] serve, in some sense, to complement this article. 
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